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ABSTRACT 

A parameterization is presented that provides extinction cross section < x e , single-scattering albedo w, and 
asymmetry parameter g of ice crystals for any combination of volume, projected area, aspect ratio, and crystal 
distortion at any wavelength in the shortwave. Similar to previous parameterizations, the scheme makes use of 
geometric optics approximations and the observation that optical properties of complex, aggregated ice 
crystals can be well approximated by those of single hexagonal crystals with varying size, aspect ratio, and 
distortion levels. In the standard geometric optics implementation used here, cr e is always twice the particle 
projected area. It is shown that w is largely determined by the newly defined absorption size parameter and the 
particle aspect ratio. These dependences are parameterized using a combination of exponential, lognormal, 
and polynomial functions. The variation of g with aspect ratio and crystal distortion is parameterized for one 
reference wavelength using a combination of several polynomials. The dependences of g on refractive index 
and ai are investigated and factors are determined to scale the parameterized g to provide values appropriate 
for other wavelengths. The parameterization scheme consists of only 88 coefficients. The scheme is tested for 
a large variety of hexagonal crystals in several wavelength bands from 0.2 to 4 /cm, revealing absolute dif- 
ferences with reference calculations of w and g that are both generally below 0.015. Over a large variety of 
cloud conditions, the resulting root-mean-squared differences with reference calculations of cloud re- 
flectance, transmittance, and absorptance are 1.4%, 1.1%, and 3.4%, respectively. Some practical applications 
of the parameterization in atmospheric models are highlighted. 


1. Introduction 

A substantial uncertainty in properly representing ice 
clouds and consequently their climate feedbacks is 
characterizing their radiative properties (Stephens et al. 
1990; Stackhouse and Stephens 1991; Schlimme et al. 
2005; Fu 2007). The fundamental radiative properties 
of atmospheric ice crystals for atmospheric models are 
the ice extinction cross section, single-scattering albedo, 
and the first moment of the scattering phase function, 
commonly referred to as the asymmetry parameter 
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(Coakley and Chylek 1975; Fu 1996; Yang et al. 2000; Fu 
2007). An increasing number of parameterizations for 
these optical properties of ice clouds are available. Such 
parameterizations generally relate the optical properties 
in selected wavelength bands in terms of predicted or 
imposed bulk characteristics of the ice, such as effective 
size, shape, and ice water content (e.g., Fu and Liou 
1993; Fu 1996; Wyser and Yang 1998; Kristjansson et al. 
1999; Y ang et al. 2000; McFarquhar et al. 2002; Key et al. 
2002; Edwards et al. 2007; Xie et al. 2012), sometimes in 
combination with atmospheric conditions, such as tem- 
perature and ice supersaturation (Kristjansson et al. 
2000; Baran 2012). Although ice crystals are found with 
a virtually countless variation of shapes (Baran 2009), 
the extinction cross section and single-scattering albedo 
of randomly oriented individual ice crystals mainly de- 
pend on their projected areas and volumes, while the 
asymmetry parameter is shown to be determined mainly 
by the aspect ratios (defined throughout as the ratio of 
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hexagonal prism height to width) of the ice crystal hex- 
agonal components and their microscale surface rough- 
ness or crystal distortion (Iaquinta et al. 1995; Macke 
et al. 1996; Fu 2007; Yang and Fu 2009; Baran 2009; van 
Diedenhoven et al. 2012a). 

Parallel to advances in understanding ice optical 
properties are efforts to improve ice microphysics mod- 
eling schemes. Some of these schemes are becoming 
increasingly complex and calculate or parameterize ice 
fall speeds and capacitances (an electrostatic analog 
used to compute ice deposition and sublimation rates) 
in terms of ice mass, projected area, and aspect ratio 
(Bohm 1989, 1992; Fleymsfield and Iaquinta 2000; 
Wood et al. 2001; Westbrook 2008; Westbrook et al. 
2008; Sulia and Harrington 2011; van Diedenhoven 
et al. 2012b). Since mass m can be related to bulk ice 
volume V via V = mlpj , where p, is the density of bulk 
ice, the ice characteristics determining fall speeds and 
capacitances are largely the same as those determining 
ice optical properties. Only crystal distortion does not 
affect fall speeds and capacitances substantially while 
being important for ice crystal optical properties. Self- 
consistency within a model dictates that the same ice 
volume, area, and aspect ratio used in an ice micro- 
physics scheme should also be used in a model’s radia- 
tive transfer scheme (Mitchell et al. 2008; Baran 2012). 
However, there is at present a lack of optical property 
parameterizations flexible enough to accommodate all 
possible combinations of ice volume, area, aspect ratio, 
and crystal distortion. Furthermore, optical properties 
parameterizations are usually given for a predetermined 
set of wavelength bands, which are not necessarily the 
same as those used in a model’s radiative transfer calcu- 
lations. Finally, most current optical property schemes 
require many and large tables to be read by the radiative 
transfer algorithm. 

In this paper, we present a simple yet flexible pa- 
rameterization that provides the extinction cross sec- 
tion, single-scattering albedo, and asymmetry parameter 
of ice crystals for any combination of volume, projected 
area, aspect ratio (of crystal components), and crystal 
distortion and at any wavelength in the shortwave. 
Similar to previous parameterizations (e.g., Fu 1996, 
2007), our scheme makes use of geometric optics ap- 
proximations and the observation that optical properties 
of complex, aggregated ice crystals can be approximated 
by those of single hexagonal crystals with varying size, 
aspect ratio, and crystal distortion. 

After describing the reference calculations, the the- 
ory, and the parameterization in section 2, we assess 
its accuracy in section 3. We provide an example of 
an application in section 4 and present conclusions in 
section 5. 


2. Theory and parameterizations 

Here, we first discuss the dependencies of extinction 
cross section, absorption cross section, and asymmetry 
parameter on ice crystal volume, projected area, aspect 
ratio, and crystal distortion before parameterizing those 
dependencies. 

a. Reference calculations 

The reference calculations are based on single hex- 
agonal plates and columns. The optical properties for 
hexagonal plates and columns with random orientation 
are calculated using the geometric optics (GO) code 
developed by Macke et al. (1996). In this standard 
geometric optics implementation, the extinction cross 
section is twice the projected area of crystals in random 
orientation. This Monte Carlo ray-tracing code takes 
crystal distortion of ice crystals into account in a statis- 
tical manner by perturbing the normal of the crystal 
surface from its nominal orientation by an angle that, for 
each interaction with a ray, is varied randomly with 
uniform distribution between 0° and 8 X 90°, where 5 is 
referred to as the distortion parameter. Thus, 8 does not 
represent a single realization of a distorted crystal but 
rather the stochastic large-scale distortion of a collection 
of ice crystals (Macke et al. 1996). For a large collection 
of ice crystals, microscale surface roughness and large- 
scale particle distortion both lead to a similar randomi- 
zation of the angles between crystal facets, which in turn 
leads to the suppression of features in the scattering 
phase matrix (Yang et al. 2008; Neshyba et al. 2013). 
Increasing the number of impurities within ice crystals 
also has a similar effect (Hess et al. 1998). Thus, the 
distortion parameter used here can be considered as 
a proxy for the randomization of the angles between 
crystal facets caused by any of these effects. 

Calculations are made for bands with the average 
wavelengths and complex refractive indices listed in 
Table 1. Various sets of calculations including different 
dimensions of the hexagonal particles are used to con- 
struct and evaluate the parameterization, as detailed in 
the next subsections. The relations between volume V, 
orientation-averaged projected area A p , maximum half- 
width of the hexagonal planes « hex . and aspect ratio a of 
hexagonal plates and columns are given in appendix A. 

The reference calculations and thus the parameteri- 
zation presented here are based on conventional geo- 
metric optics calculations. It is generally accepted that 
GO approximations are applicable to scattering size 
parameters above about 100 (e.g., Macke et al. 1995; 
Yang et al. 2000, 2004, 2013). (Here we define the particle 
scattering size parameter y- sca , as 2-nT sph /A, where r sph 
is the radius of a sphere with the equivalent projected 
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Table 1. Average wavelengths A„ and real ( m r ) and imaginary 
(m,-) parts of the refractive index used in the reference calculations. 
Values of m n m , (Warren 1984; Warren and Brandt 2008), and A„ are 
weighted by solar irradiance at the top of the atmosphere within each 
wavelength band. Band edges are given in Table 1 of Kato et al. 
(1999), with bands 9-15 there corresponding to band 9 here. 

A„ (/im) 

m r 

rrii 

0.256 

1.3480 

8.082 X 10~ 9 

0.280 

1.3407 

6.751 X 10~ 9 

0.296 

1.3353 

5.756 X 10 -9 

0.319 

1.3307 

4.878 X 10~ 9 

0.335 

1.3275 

4.269 X 10~ 9 

0.365 

1.3231 

3.420 X 10~ 9 

0.420 

1.3177 

2.261 X 10~ 9 

0.482 

1.3140 

1.742 X 10~ 9 

0.598 

1.3098 

7.511 X 10~ 9 

0.690 

1.3071 

2.445 X 10~ 8 

0.719 

1.3056 

7.549 X 10“ 8 

0.762 

1.3065 

3.834 X 10~ 8 

0.813 

1.3047 

1.376 X 10“ 7 

0.862 

1.3038 

2.330 X 10“ 7 

0.926 

1.3028 

5.267 X 10“ 7 

1.005 

1.3014 

1.724 X 10~ 6 

1.111 

1.2997 

2.228 X 10 -6 

1.333 

1.2955 

7.335 X 10“ 5 

1.562 

1.2906 

4.841 X 10” 4 

1.770 

1.2837 

2.627 X 10” 4 

2.051 

1.2717 

1.212 X 10“ 3 

2.210 

1.2629 

3.064 X 10" 4 

2.584 

1.1815 

2.773 X 10“ 2 

3.284 

1.4310 

2.719 X 10 _1 

3.809 

1.3874 

7.558 X 10~ 3 

4.292 

1.3473 

1.639 X 10~ 2 


area; i.e., r sp h = \jA p hr.) For smaller sizes, the errors on 
ice crystal optical properties and ice cloud radiation cal- 
culations caused by GO approximations are not yet ad- 
equately quantified. Such a quantification would require 
intensive calculations using techniques that are able to 
accurately calculate optical properties of a variety of ice 
crystal shapes with small size parameters (e.g., Yang et al. 
2013). Such an investigation is beyond the scope of the 
current paper and left for future work. 

b. Extinction cross section and single-scattering 
albedo 

The single-scattering albedo of a single ice particle is 
defined to be 


(0 = 1 -°^, ( 1 ) 

a e 

where the extinction cross section of a particle is given by 

<re = QAp’ ( 2 ) 


and Liou 1996). According to anomalous diffraction 
theory (ADT; van de Hulst 1957; Bryant and Latimer 
1969; Mitchell and Arnott 1994), the absorption cross 
section can be approximated by 


1 — exp 


—A'rrm i V 
A A.. 


A r ’ 


(3) 


where m, is the imaginary part of the refractive index 
(Table 1) at wavelength A, and V is the volume of the ice 
crystal, which is related to mass m via V = ml pi. From Eqs. 
(l)-(3), it is apparent that the single-scattering albedo at 
a given wavelength is mainly determined by the ice crystal 
effective distance or photon path d e as given by Bryant and 
Latimer (1969): 



The relation between the single-scattering albedo of a 
compact hexagonal particle with aspect ratio a of unity 
and the effective distance for the wavelength bands with 
substantial absorption is shown in Fig. la. For this, calcu- 
lations were made for smooth crystals with 100 different 
ratios of volume to projected area varying between 1 
and 135 pum. The effect of particle distortion on single- 
scattering albedo is considered negligible (Yi et al. 2013). 

The ice crystal effective distance is related to the 
commonly used, but variably defined (McFarquhar and 
Heymsfield 1998) bulk ice effective diameter D e , in 
which “bulk” denotes an average over a particle size 
distribution. Flere we define D e through the ratio of the 
size-distribution integrated volume and project area: 


D 


e 


3 

2 


*00 

Jo 


V(D)N(D) dD 


A (D)N(D) dD 
o p 


(5) 


where N is the number density of the ice crystals of 
maximum dimension D (Foot 1988; Francis et al. 1994; 
Baum et al. 2005a). The smooth relation of the single- 
scattering albedo to the ratio of volume to area has 
often been used to parameterize the single-scattering 
albedo as a function of bulk ice effective diameter (e.g., 
Wyser and Yang 1998; Yang et al. 2000; McFarquhar 
et al. 2002; Key et al. 2002; Edwards et al. 2007; Xie et al. 
2012 ). 

However, Eq. (3) suggests that a more fundamental 
parameter governing the single-scattering albedo of in- 
dividual particles is 


where Q e is the extinction efficiency, which is equal to 2 
in the geometric optics limit (Macke et al. 1996; Yang 


Labs 


m -V m i 

lA- p = T d *' 


( 6 ) 
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Fig. 1. Single-scattering albedo for compact hexagonal ice crystals as a function of (a) effective distance and 
(b) absorption size parameter. Colors indicate different wavelength bands. The gray-red line represents the expo- 
nential fit described in the text. The blue dashed line is the result using ADT. 


which we will refer to as the absorption size parameter. 
Note that ^ abs is dimensionless. As seen in Fig. lb, the 
single-scattering albedo in nearly all bands collapse to 
a single function of ;y abs . Only single-scattering albedos 
corresponding to the two bands with the strongest ab- 
sorption (2.6 and 3.3 gm with m, > 0.02 as seen in Table 1) 
substantially deviate from the single curve. For the re- 
maining bands, single-scattering albedos at the same ab- 
sorption size parameter but for different wavelengths 
agree within about 0.01, which is discussed in greater 
detail later. The dependence of single-scattering albedo 
for compact particles (o> a=1 ) to absorption size parameter 
is found to be fit well by an exponential function: 

" a =i = 1 - « n 0 ~ < r<w,bs )> ( 7 ) 

where the form of the parameterization guarantees that 
oj,, = ] = 1 for , y abs = 0. Using a least squares fit, we find 
a 0 = 0.457 593 and a i = 20.9738. (The parameterization 
coefficients are computed in single precision with 4-byte 
floating points and given to six significant digits through- 
out.) The values corresponding to the 2.6- and 3.3-/xm 
bands with strong absorption were excluded when per- 
forming this fit. The root-mean-square difference be- 
tween the fit and the single-scattering albedos used for 
the fit is 0.002. As evident from Fig. lb, errors for strongly 
absorbing wavelengths are larger, as further discussed in 
section 3. As can be derived from Eqs. (l)-(3), in con- 
ventional ADT, the coefficients ci 0 and cq of Eq. (7) are 
l/Q e and 4 t t, respectively (Mitchell and Arnott 1994). 
However, Fig. lb shows that ADT generally over- 
estimates the single-scattering albedo of compact hex- 
agonal crystals. As discussed by Mitchell et al. (1996) 
among others, this error is principally attributable to the 


neglect of internal reflections. Several modifications to 
ADT have been proposed (e.g., Mitchell 2002; Yang 
et al. 2004), which generally include several tuning 
factors; an evaluation of such schemes is beyond the 
scope of this study. Note that a quantity related to the 
absorption size parameter [Eq. (6)] was used previ- 
ously by Takano and Liou (1989) to parameterize co for 
weakly absorbing (co < 0.75) hexagonal crystals at 
wavelengths shorter than 2.5 fjem. 

The parameterization described by Eq. (7) assumes 
compact hexagonal ice crystals ( a = 1). When aspect 
ratios deviate from unity by about a factor of 2, co is 
generally greater than that of a compact hexagonal 
crystal, as seen in Fig. 2. This effect is stronger for plates 
(Fig. 2a) than for columns (Fig. 2b). Presumably this 
increase of co with aspect ratio deviating from unity re- 
sults from a decrease of the orientation-averaged path- 
length through the particle, because internal reflections 
are maximized for compact particles. Similar to Fig. lb, 
Fig. 2 shows that, for a given aspect ratio, co collapses to 
a single function of Note however that slight dif- 
ferences in co at the same ^- abs but for different combi- 
nations of d e and wavelengths result in “spikes” in Fig. 2. 
Such spikes may result from wave resonance effects, 
leading to rays orbiting the particle through a series of 
internal reflections for specific combinations of wave- 
length, refractive index, and particle size, thereby in- 
creasing absorption (Guimaraes and Nussenzveig 1992; 
Mitchell et al. 2006). Interestingly, the spikes diminish as 
aspect ratio increasingly deviates from unity, which is 
consistent with such ray-orbiting phenomena (Mitchell 
et al. 2006). However, some Monte Carlo noise is also 
apparent in Fig. 2, which increases with aspect ratio 
deviation from unity. Figure 2 also shows functional fits 
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Xobs 


FIG. 2. As in Fig. lb, but for (a) hexagonal plates and (b) columns 
with various aspect ratios, as indicated by the colors. All bands 
listed in Table 1 with weak to moderate absorption (m ; < 0.02) are 
included. Gray-red dashed lines indicate fits. 



Xobs 


Fig. 3. Differences between the parameterized single-scattering 
albedo for particles with an aspect ratio of unity and the calculated 
single-scattering albedos for several (a) plates and (b) columns with 
various aspect ratios. Gray-red dashed lines indicate fitted log- 
normal functions as described in the text. 


to the ^ a bs versus- w relations for plates with aspect ratios 
of 0.8, 0.5, 0.2, 0.1, 0.05, and 0.02 and for columns with 
aspect ratios that are the inverse of those for the plates. 
First, the difference Aw„ between the calculated single- 
scattering albedo and the exponential fit assuming 
compact hexagonal crystals [Fig. lb and Eq. (7)] is 
determined for each of the aspect ratios, as shown in 
Fig. 3. We find that A co a can be approximated by log- 
normal functions that are also shown in Fig. 3; that is, 


Aw 


\J2ttI 


1 Tibs 


-exp 


0"*abs-/2) 


2/f 


( 8 ) 


The fit parameters /,-, where i = 0-2, are determined for 
each of the aspect ratios using a least squares fitting 
method. In turn, /, are found to be smooth functions of 
the logarithm of a that, for either plates and columns, 
can be approximated by cubic functions; that is. 


*,■(«)= Zc.Tog >a, (9) 

;'=o 

for which the coefficients Cy are listed in Table 2. 
(Throughout, “log” refers to base 10.) To better con- 
strain this latter fit, we also include calculations for as- 
pect ratios of 0.01 and 100, which are not shown in Figs. 2 
and 3. Finally, the single-scattering albedo of a hexago- 
nal ice crystal of any size and any aspect ratio can be 
estimated for any wavelength by 

w CW“) = £t W +Aa V ( 10 ) 

where u > a= i and A co a are calculated using Eqs. (7) and 
(8), respectively. Note that the spikes and noise seen 
in Fig. 2, as mentioned above, are even more apparent in 
Fig. 3. Given that such spikes are about 0.01 u> units in 
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Table 2. Coefficients C/j used by Eq. (9) to compute coefficients /,- 


in Eq. (8). See also Fig. 4, box 3. 


j 

C 0j 

c u 

C 2j 

Plates 

0 

0.000527060 

0.309748 

-2.58028 


1 

0.008675 96 

-0.650188 

-1.34949 


2 

0.0382627 

-0.198214 

-0.674495 


3 

0.010855 8 

-0.0356019 

-0.141318 

Columns 

0 

0.000378774 

0.390452 

-2.36821 


1 

0.004632 83 

0.420040 

1.07603 


2 

0.00593106 

-0.0848059 

-0.729 980 


3 

-0.001 17167 

0.018 6601 

0.232446 


magnitude, the parameterization of co (see Fig. 4 for 
a summary) is limited to an overall accuracy of about 
0.01 at wavelengths where m t < 0.02. The resulting pa- 
rameterized single-scattering albedo as a function of the 
absorption size parameter and aspect ratios for plates 
and columns is shown in Fig. 5. 

As given by Eq. (2), the extinction cross section is 
always twice the projected area of crystals in random 
orientation in the standard geometric optics implemen- 
tation used here (see Fig. 4, box 5). This approximation 
generally holds for scattering size parameters above 
about 50 (e.g., Yang et al. 2000, 2004, 2013). For particles 
with smaller scattering size parameters (1-50), so-called 
edge effects and wave interference have to be taken into 
account (Mitchell 2000; Yang et al. 2013), which lead to 
a general increase of the extinction coefficient up to about 


3 with decreasing size parameter and to oscillations in 
the extinction coefficient in size parameter space of about 
20% or smaller. For yet smaller size parameters (;^ scat < 1), 
particles are comparable to or smaller than the wave- 
length and transition to the Rayleigh regime, associated 
with a decrease of extinction coefficient with decreasing 
size. Since the parameterizations presented in this paper 
are based on calculations that are based on the standard 
geometric optics approximations, these effects are ne- 
glected and Q e is assumed to be equal to 2 for all sizes. 
Thus, substantial errors in the extinction coefficient ow- 
ing to the geometric optics approximation can be ex- 
pected for small ice crystals. However, the contribution of 
these small crystals to the total extinction cross section is 
expected to be generally minor under many conditions, 
owing to their small projected area. To our best knowl- 
edge, a proper quantilication of the errors in extinction 
cross sections and ice cloud radiation simulations result- 
ing from the approximation Q e = 2 is not yet available. 
As discussed in section 2a, such an investigation is beyond 
the scope of the current paper and left for future work. 

c. Asymmetry parameter 

In the geometries optics regime, the total scattering 
phase function Plot includes contributions from internal 
refraction and reflections P RT (determined using ray 
tracing) and contributions from diffraction /' d ii (Macke 
et al. 1996): 



Fig. 4. Summary of the parameterization of single-scattering albedo (box 4), with the absorption size parameter in 
box 1, parameterization of single-scattering albedo for particles with aspect ratio of unity (&t 0 =i) in box 2, and the 
parameterization of the relative difference between cu a=1 and scattering albedo for particles with other aspect ratios 
in box 3. The geometric optics approximation for the extinction cross section is given in box 5. 
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Cl) 


Plates 0.50 0.58 0.67 0.75 0.83 0.92 1.00 Columns 



0.1 1.0 10 
Aspect ratio 

Fig. 5. Parameterized single-scattering albedo as a function of particle aspect ratio and absorption size parameter. 


P UA (0)=^[(2a>-l)P RT (0) + P' m m- 


we simply approximate it with a parametric form justi- 
fled by the integral over a diffraction pattern: 


The scattering phase function is often characterized by 
its first moment, the asymmetry parameter g tot (Macke 
et al. 1996): 


Slot 


r i r 

P t . cos 9 sinfl d6 = — [(2 co 
, o 2 co 


1)#RT + £dif] ■ 

(12) 


where g dif and g RT are the asymmetry parameters of the 
diffractive contributions and the refractive plus reflec- 
tive (or ray tracing) contributions to the phase function, 
respectively. A parameterization of co is presented in 
the previous section. Note that at nonabsorbing wave- 
lengths ( co = 1), Eq. (12) reduces to 

Stot = ^(gRT + Sdif)- ( 13 ) 


Here we aim to parameterize g dii and g RT as a function 
of particle volume, projected area, aspect ratio, distor- 
tion parameter, and wavelength. Then Eq. (12) can be 
used to obtain the total asymmetry parameter for any 
hexagonal ice crystal at any wavelength. 

As shown in Fig. 6, the diffraction asymmetry pa- 
rameter is mainly determined by the particle scattering 
size parameter ;y scat . For ;y scat larger than about 100, 
gdif = 1- For smaller ^ SC at, gdii can be calculated analyt- 
ically with series of Bessel functions (van de Hulst 1957; 
Macke et al. 1996; Mishchenko and Macke 1998). Here, 


gdif = V ftlln * s “‘ + b 2> ( 14 ) 

with b 0 = -0.822 315, b t = -1.201 25, and b 2 = 0.996 653, 
as shown in Fig. 6. The coefficients were obtained using 
a least squares fit. The diffraction asymmetry parameter 



Fig. 6. Diffraction asymmetry parameter of hexagonal ice crys- 
tals with an aspect ratio of unity as a function of scattering size 
parameter. The gray dashed line represents the exponential fit 
described in the text. 
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Fig. 7. Summary of the parameterization of asymmetry parameter g t ot (box 9), with parameterization of the diffraction asymmetry 
parameter in boxes 1 and 2, parameterization of ray-tracing asymmetry parameter at 862 nm in boxes 3-5, and correction factors for real 
refractive index and absorption in boxes 6-8. 


is limited to be greater than 0.5. This parameterization of 
gun is summarized in Fig. 7, boxes 1 and 2. 

The refraction-plus-reflection asymmetry parameter 
depends on particle volume, projected area, aspect ratio, 
distortion parameter, and complex refractive index. 
However, at any given nonabsorbing wavelength, g RT of 
single hexagonal ice crystals mainly depends on aspect 
ratio and crystal distortion (cf. Fu 2007; Yang et al. 2008; 
Yang and Fu 2009; van Diedenhoven et al. 2012a). For 
example. Fig. 8 shows the total asymmetry parameter at 
a nonabsorbing wavelength of A = 862 nm as a function 
of 8 for several aspect ratios. Note that these calculations 
are made for large particles (^ scat 100), and thus g dil = 1 
and the dependencies of the total asymmetry parameter 
on aspect ratio and distortion parameter seen in Fig. 8 
are solely attributable to variations in g RT [see Eq. (13)]. 
Here, the wavelength band at 862 nm is used as a refer- 
ence because it corresponds to similar bands commonly 
used in aircraft and satellite retrievals of ice cloud 
properties (e.g., Chepfer et al. 2001; King et al. 2004; 
van Diedenhoven et al. 2012a, b). As described by van 
Diedenhoven et al. (2012a), sizes are varied so that the 
projected areas of the particles, assuming random ori- 
entation, correspond to the projected areas of spheres 
with radii of 14, 20, 28, 40, 56, 80, 113, 160, 226, and 
320 pm. The aspect ratio of columns is varied between 
1 and 50 with 26 geometrically increasing steps. The as- 
pect ratios of plates are the inverse of those for columns, 
for a total of 51 aspect ratios. The distortion parameter is 


varied between 0 and 0.8 in steps of 0.05. Since GO 
calculations are independent of size for large particles, 
results are simply averaged over all sizes in order to 
decrease numerical noise. In the following, we aim to 
first parameterize the total asymmetry parameter at 
A = 862 nm from which we derive g RT at A = 862 nm 
using Eq. (13). Subsequently, we derive factors to scale 



Fig. 8. Ice crystal total asymmetry parameters at 862 nm as 
a function of distortion parameter for various aspect ratios. Dashed 
lines represent the fits described in the text. 
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the resulting g RT to other wavelengths — that is, for other 
complex refractive indices. 

Figure 8 shows that asymmetry parameters increase as 
aspect ratio departs from unity owing to the increase 
of parallel surface areas, leading to greater probability 
of light passing through the particle with low orders of 
refraction plus reflection and a minimal change of di- 
rection (Yang and Fu 2009). Increase of crystal distor- 
tion increases the chance of light refracted away from 
forward directions, leading to decreased asymmetry 
parameters (Yang et al. 2008). 

Figure 8 shows that the dependence of g a=1 on 5 
for compact particles (a = 1) can be approximated by 
a fourth-degree polynomial in S with coefficients p 0 = 
0.780550, Pl = 0.00510997, p 2 = -0.087 8268, p 3 = 
0.111 549, and p 4 = —0.282453. Similar to the procedure 
followed above in the parameterization of single- 
scattering albedo, the differences A g a between g a=l 
and g tot at all other aspect ratios in our database are then 
fit with third-degree polynomials in 8 with coefficients 
p A In turn, the dependence of those p A i values on login- ) 
is found to be well approximated by sixth-degree poly- 
nomial functions in a for both columns and plates. (The 
degrees of the polynomials are chosen to keep the max- 
imum differences of parameterized and reference g tot 
values below 0.005.) The parameters for this fit are listed 
in Table 3, and the results are shown in Fig. 8. Thus, A g a is 
determined by first computing from a the p A i values, and 
subsequently computing from 5 the A g„ values for the 
required a-8 combination. In turn, the total asymmetry 
parameter at 862 nm for the required values of aspect 
ratio and distortion is computed from 

g tot (a,5, A = 862 nm) =g a=1 (S) + Ag a (a,5). (15) 

This parameterization of the total asymmetry parameter 
at A = 862 nm is summarized in Fig. 7, boxes 3-5. 

Figure 9 shows the parameterized g tot at 862 nm and 
the differences from the reference calculations (cf. Fig. 1 
in van Diedenhoven et al. 2012a). The asymmetry pa- 
rameter at 862 nm is seen to be approximated with 
an accuracy better than 0.004. The refraction plus re- 
flection asymmetry parameter at 862 nm can be derived 
from this parameterization using g RT = 2g tot — 1, since 
gdif = 1 and co = 1 [see Eq. (13)]. Next, we will dem- 
onstrate how this parameterized asymmetry parameter 
for A = 862 nm can be scaled to other wavelengths — that 
is, for other complex refractive indices. 

The dependence of g RT on the real part of the re- 
fractive index m r is investigated by performing geo- 
metric optics calculations for smooth large hexagonal 
crystals at the wavelength bands listed in Table 1, but 
artificially setting the imaginary part of the refractive 


Table 3. Coefficients of the sixth-order polynomial to compute 
coefficients p A l in Eq. (15). See also Fig. 7, box 4. 



i 

loj 

J 

<12J 

Plates 

0 

-0.001 33106 

-0.000782076 

0.002 05422 


l 

0.0408343 

-0.00162734 

0.024092 7 


2 

0.525289 

0.418336 

-0.818352 


3 

0.443 151 

1.53726 

-2.403 99 


4 

0.008 525 15 

1.88625 

-2.64651 


5 

-0.123 100 

0.983 854 

-1.291 88 


6 

-0.037 6917 

0.187708 

-0.235 359 

Columns 

0 

-0.001 89096 

0.000637 430 

0.001 573 83 


1 

0.00981029 

0.0409220 

0.00908004 


2 

0.732 647 

0.053 9796 

-0.665 773 


3 

-1.59927 

-0.500870 

1.86375 


4 

1.54047 

0.692547 

-2.053 90 


5 

-0.707187 

-0.374 173 

1.012 87 


6 

0.125 276 

0.0721572 

-0.186466 

index to 

zero 

(i.e., co = 1). 

Figure 10 shows that g RT 

relative to its value at 862 nm (m r = 1.3038) is a smooth 

function of m r 

for a given aspect ratio. We hypothesize 


that, as Fresnel reflection on the particle surfaces is in- 
creased with an increasing refractive index, forward 
scattering is reduced and backward scattering is in- 
creased (i.e., asymmetry parameter is lowered). It is seen 
in Fig. 10 that the ratio of the asymmetry parameter with 
a certain refractive index m r ] to that with another re- 
fractive index m r2 is reasonably well approximated by 
functions of the form 

_K i -e)(m r 2 + e) 
m ' (m rl +e)(m r2 -ey 1 ' 

where e is a fitting coefficient that depends on aspect 
ratio. Since we use calculations at 862 nm as a reference — 
here m rl = 1.3038. The functional form of Eq. (16) is 
inspired by the Fresnel equation for reflection (Born and 
Wolf 1999). The parameter e is computed from the value 
that minimizes the root-mean-square difference between 
Eq. (16) and gwripirg) over the range of m r2 in Table 1. 
The value of e is found to be 0.960251 for compact ice 
crystals (a = 1), and it decreases as the aspect ratio in- 
creasingly deviates from unity. We find that e scales well 
with log(o-); that is, 

e(«) = c 0 + Ci logo: , (17) 

with c 0 = 0.960251 and c 4 = 0.429 181 for plates and c 0 = 
0.941791 and c\ = —0.216010 for columns. Essentially 
the same fit parameters are found when taking other m r l 
values and for particles with increased 8 (not shown). 
This parameterization of dependence of g RT on m r is 
summarized in Fig. 7, box 6. 


1772 


JOURNAL OF THE ATMOSPHERIC SCIENCES 


Volume 71 


Asymmetry parameter 



0.1 1.0 10 
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Fig. 9. (top) Parameterized total asymmetry parameters at 862 nm and (bottom) the differences from the reference 

calculations. 


As noted by Fu (2007), Eq. (12) shows that the de- 
pendence of the total asymmetry parameter on single- 
scattering albedo is largely determined by the weighting 
of the ray-tracing and diffractive parts of the asymmetry 
parameter. However, g RT itself also depends on the 
single-scattering albedo; ignoring this dependence can 
lead to errors in the total asymmetry parameter up to 
about 5%, depending on the wavelength band and as- 
pect ratio, as shown in Fig. 11a. This effect was neglected 
by Fu (2007). To investigate the dependence of g RT on 


absorption, Fig. lib shows g RT relative to its value when 
co = 1 as a function of coalbedo (1 — co). The value of 
g RT is found to increase nearly linearly with (1 — co) for 
(1 — co) s 0.3, but to decrease for greater (1 — co). This 
dependence of g RT on co is seen to vary with aspect ratio. 
The dependence also varies with the real part of the 
refractive index, as evident from the differences be- 
tween wavelengths above and below 3/rm. The de- 
pendence is not found to be substantially affected by 
increasing distortion (not shown). We hypothesize that 
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Fig. 10. The refraction-plus-reflection asymmetry parameter of 
nonabsorbing ice crystals as a function of real refractive index, 
relative to its value at 862 nm, or m r = 1.3038. Calculations for 
several plate aspect ratios are indicated by color. Dashed lines 
represent fits as described in the text with indicated fit parameters e. 

the dependence of g RT on w results from two competing 
effects with increasing absorption: 1) contributions of 
shorter paths through the particle with lower orders of 
scattering, which generally are directed forward (Yang 
and Fu 2009), are increased, tending to increase g RT , and 
2) contributions of rays reflected sideways or backward 
on the particle surface (i.e., external reflections) are in- 
creased relative to the rays entering and exiting the 
particle, tending to decrease g RT . To determine a factor 
to partly correct the refraction plus reflection asymmetry 
parameter that was parameterized above assuming no 
absorption [i.e., g RT (o> = 1)] for variations in coalbedo, 
we first fit a fifth-degree polynomial in (to — 1) through 
the calculated g R - r /g R r (oj = 1) values for a = 1 shown in 
Fig. lib. For simplicity, here we do not account for the 
dependence of g RT on real refractive index. The six 
coefficients s t for this fitted correction factor C w> i are 
1.00014, 0.666094, -0.535 922, -11.7454, 72.3600, and 
— 109.940. Next, the relative differences (denoted 
between C,„j and calculated values of gRx/girK" = 1) 
for aspect ratios deviating from unity are determined 
and fit by the linear function C w 2 = k(a> — 1) + 1, where 
A: is a fit parameter that depends on aspect ratio. Analo- 
gous to the procedure used above, we finally fit the de- 
pendence of k on aspect ratio using the linear functional 
form k = it log(cr). Thus, C„, i2 is computed by 

C a)2 = u log(a)(*;-l) + l, (18) 

where u = —0.213 038 for plates and u = 0.204016 for 
columns. For illustration, the fit result for a = 0.1 is also 
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Fig. 11. (a) Relative differences between calculated total 
asymmetry parameters and those estimated by using Eq. (12) and 
the refraction-plus-reflection asymmetry parameter for non- 
absorbing particles (to = 1). Calculations for different wavelength 
bands are indicated by colors. Filled circles represent compact 
particles (a = 1) and open diamonds indicate thin plates (a = 0.1). 
Only two aspect ratios are shown to reduce clutter, (b) Relative 
differences in the refraction-plus-reflection asymmetry parameter 
for particles with w = 1 and values for particles with lower single- 
scattering albedos. Red lines indicate fits for a = 1 (dashed) and 
a = 0.1 (dashed-dotted) as described in the text, (c) As in (a), but 
using the derived factors as described in the text to correct the 
refraction-plus-reflection asymmetry parameter for absorption 
effects. 
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0.5 0.6 0.7 0.8 0.9 1.0 

£j calculated 


FIG. 12. Absolute difference between parameterized and calcu- 
lated single-scattering albedos for several wavelength bands as 
indicated by colors. Circles, diamonds, and squares indicate hex- 
agonal columns, plates, and compact hexagons with aspect ratio of 
unity, respectively. 


shown in Fig. lib. The computation of these factors to 
partly correct g RT (w = 1) for variations in coalbedo is 
summarized in Fig. 7, boxes 7 and 8. Figure 11c shows 
the difference between the calculated g tol values and 
those parameterized using g R x(w = 1) with correction 
factors C a) ] and C„ )2 - The errors are generally below 1% 
for A < 2.2 /rm and below 2% for longer wavelengths. 

Finally, the asymmetry parameter for hexagonal crys- 
tals can be approximated via 

g tot (a,8, m r ,u)) 

1 

= Ico ~ + £difJ ’ 

(19) 

as summarized in Fig. 7, box 9. The single-scattering 
albedo a> is parameterized as described in section 2b 
and Fig. 4. 

3. Results 

a. Optical properties of individual particles 

The parameterization is tested on GO calculations for 
hexagonal crystals with aspect ratios varying between 
0.05, 0.1, 0.5, 1, 2, 5, and 20; effective distances of 10, 30, 
50, 70, 90, 110, 150, and 190 gm; and distortion parameters 
of 0, 0.2, 0.4, and 0.6 and at wavelength bands listed in 
Table 1. For all wavelength bands with substantial ab- 
sorption, absolute differences between calculated and pa- 
rameterized single-scattering albedos are generally below 
0.015, as seen in Fig. 12. Only for bands with m, > 0.02 



0.70 0.75 0.80 0.85 0.90 0.95 1.00 

g calculated 

Fig. 13. As in Fig. 12, but for asymmetry parameter. 


(at 2.6 and 3.3 gm) are absolute differences up to 0.05 
obtained, especially for small single-scattering albedos, 
which is expected considering the different dependence 
of co on ^ abs for these bands seen in Fig. 1. Though not 
evident in Fig. 12, the smallest errors are obtained for 
relatively compact particles with 0.5 < a < 2. 

Figure 13 shows the absolute differences between 
calculated and parameterized asymmetry parameters 
for all wavelength bands listed in Table 1. For non- 
absorbing bands, the accuracy of the parameterization is 
better than 0.01. For moderately absorbing wavelengths, 
absolute differences between calculated and parame- 
terized asymmetry parameters do not exceed 0.015, but 
the parameterization is slightly biased toward under- 
estimating the asymmetry parameter, especially for 
larger, more absorbing particles. For strongly absorbing 
wavelengths (m ; > 0.02), errors in the parameterized 
asymmetry parameter up to 0.05 are seen, which is at- 
tributable to the deviations of single-scattering albedo 
from the parameterization for these bands, and the fact 
that the correction of g RT for its dependence on the 
single-scattering albedo for these bands are not well 
represented by the fits (Fig. 11). The smallest errors 
are again obtained for relatively compact particles with 
0.5 < a < 2. 

b. Application to flux calculations 

To assess the errors on flux in heating calculations 
caused by errors in parameterized optical properties, we 
apply two-stream calculations (see appendix B; Coakley 
and Chylek 1975; Wiscombe and Grams 1976) to uni- 
form ice cloud layers. We use the same optical proper- 
ties of single ice particles throughout. To obtain optical 
properties for realistic ice size distributions, gamma size 
distributions of the form 
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Fig. 14. Relative parameterization errors in transmittance, re- 
flectance, and absorptance of uniform cloud layers. The box-and- 
whisker plots show median values, 25th and 75th percentiles, and 
full ranges of the errors in several transmittance, reflectance, and 
absorptance bins. Cloud optical thicknesses range from 0.01 to 100 
and wide ranges of different size distributions, aspect ratios, and 
distortion parameters are included, as described in the text. A solar 
zenith angle of 60° is assumed. Further details provided in text. 

^L{D) = N 0 D^e- x r D (20) 

are applied (Mitchell 1991; Heymsfield et al. 2002; Baum 
et al. 2005b), where D is the maximum particle dimension 
and ijl p and \ p are the shape and slope parameters of the 
size distribution. Here N 0 is a normalization factor. A 


range of size distributions is applied by varying \ p in 
Eq. (20) from 20 to 200 cm - 1 with increments of 20 cm~ 1 
and using the relation between /jl p and \ p determined by 
Heymsfield et al. (2002) from in situ measured ice size 
distributions; namely, fi p = A°' 64 — 2 (with \ p in units of 
per centimeter). [We note that there is ample uncertainty 
in this relation because of unknown shattering issues in 
these measurements (Korolev et al. 2011), but our results 
are largely insensitive to this assumption.] Two-stream 
calculations are made for 10 different optical thicknesses 
ranging from 0.01 to 100 and a solar zenith angle of 60°. 
Spectral reflectance and transmittance calculations are 
made for the wavelength bands listed in Table 1. Relative 
solar reflectance R and transmittance T are obtained 
through integration of the spectral quantities, weighted 
by the relative contribution to incoming spectral solar 
flux by each wavelength band (cf. Fu 1996). Absorptance 
is then defined as 1 — R — T. 

Figure 14 shows errors in reflectance, transmittance, 
and absorptance resulting from the parameterization. 
Relative errors in reflectance generally increase with 
decreasing reflectance and the same is generally true for 
the error in transmittance. The error in absorptance 
peaks around absorptance values of 0.1, which occur 
mostly at optical thicknesses of 20-60. Errors in re- 
flectances are mostly caused by errors in the parame- 
terization of asymmetry parameter, while errors in 
transmittance and absorptance are mostly affected by 
parameterization errors in the single-scattering albedo, 
as might be expected. Overall median errors and inter- 
quartile ranges in reflectance, transmittance, and ab- 
sorptance are listed in Table 4 for several aspect ratio 
ranges. Lowest errors and interquartile ranges are gen- 
erally obtained for relatively compact crystals (0.5 < 
a < 2), while more extreme aspect ratios yield somewhat 
larger errors. The overall relative root-mean-squared 
errors in reflectance, transmittance, and absorptance are 
1.4%, 1.1%, and 3.4%, respectively. Crystal distortion is 
found to have no significant effect on the resulting errors 
(not shown). Similar errors are also obtained for other 
solar zenith angles (not shown; cf. Fu 1996). 

To put our results in perspective with those of other 
parameterizations, we note that the averaged absolute 
values of the relative errors in reflectance, transmittance, 
and absorptance of 1.1%, 0.7%, and 2.5%, respectively, 


Table 4. Medians and interquartile ranges (within parentheses) of relative errors in fluxes and absorptance owing to the proposed 
parameterization for several aspect ratio ranges. All values are percentages. 


Quantity 


All a 

0.5 < a rs 2 

a < 0.1 

a > 10 

Reflectance 


0.69 (1.1) 

0.50 (0.52) 

1.4 (2.1) 

0.85 (1.0) 

Transmittance 


0.00 (0.74) 

0.00 (0.49) 

-0.07 (1.26) 

0.09 (1.15) 

Absorptance 


-1.9 (2.2) 

-1.4 (0.61) 

-2.7 (1.7) 

-4.3 (3.3) 
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resulting from our parameterization, are comparable to 
the corresponding values obtained by the parameteri- 
zation of Fu (1996), which are 1.2%, 0.3%, and 2.9%, 
respectively, for a solar zenith angle of 60° (their Fig. 6). 
The absolute maximum errors in reflectance, trans- 
mittance, and absorptance of 6%, 4%, and 11%, re- 
spectively, are also comparable or somewhat larger 
than the corresponding maximum errors reported by 
Fu (1996), which are about 6%, 2%, and 8%. Flowever, 
the parameterization scheme presented here is evaluated 
over a much wider range of particle size distributions, 
cloud optical thicknesses, aspect ratios, and distortion 
values than the scheme by Fu (1996). Moreover, our 
parameterization scheme is much more flexible than 
previous schemes. 

Finally, we note that the errors given here do not in- 
clude any possible errors owing to the application of the 
standard geometric optics approximations — the most 
important being the approximation Q e = 2. As stated 
above, quantification of such errors is currently lacking 
and left for future work. 

4. Practical application in atmospheric models 

It can now be considered conventional wisdom that 
shortwave optical properties of complex aggregates 
of hexagonal ice crystals at a given wavelength are 
mainly determined by their effective distance, distor- 
tion, and the aspect ratios of their components (e.g., 
Fu 1996, 2007; Um and McFarquhar 2007, 2009; van 
Diedenhoven et al. 2012a). Although the parameteri- 
zation proposed here is based on single hexagonal 
plates and columns, it can provide optical properties of 
particles with any combination of volume, projected 
area, aspect ratio, and distortion parameter. In the case 
of complex aggregates, such as bullet rosettes and ag- 
gregates of columns or plates, the aspect ratio to be 
used in the parameterization is the mean aspect ratio of 
the components of the complex particle, as described 
by Fu (2007). For example, in the case of bullet ro- 
settes, which are commonly defined as a structure of 
a number of bullets (usually six) that all have the same 
aspect ratio (Iaquinta et al. 1995; Yang and Liou 1998; 
Garrett et al. 2001), the aspect ratio to be used in the 
parameterization is the aspect ratio of the individual 
bullets, while the volume and projected area are quanti- 
ties determined for the whole rosette. These assumptions 
are essentially the same as those made by the parame- 
terization by Fu (1996, 2007). 

The flexibility of our parameterization can be exploi- 
ted to efficiently obtain ice optical properties that are 
consistent with microphysical assumptions made in a 
cloud or climate model. Models of cloud microphysics 


commonly use power laws to relate ice crystal mass m to 
crystal maximum dimension D: 

m = a rn Dbm ■ ( 21 ) 

Flabit-dependent power laws relating particle projected 
area to maximum dimension are also available in the 
literature (e.g., Mitchell et al. 1996; Baker and Lawson 
2006): 

A p =a A D bA (22) 


Furthermore power-law relationships are available that 
relate particle aspect ratio to maximum dimensions 
(e.g., Auer and Veal 1970; Um and McFarquhar 2007, 
2009) that can be used in treatments of cloud micro- 
physics that require assumptions about particle aspect 
ratio (e.g., Bohm 1989, 1992; Wood et al. 2001; Westbrook 
et al. 2008; Westbrook 2008). 

The parameterization presented here can be used to 
compute the ice crystal shortwave optical properties that 
are consistent with the assumed (or predicted) particle 
mass, area, and aspect ratios. Crystal distortion is the 
only parameter important for ice shortwave optical 
properties that is not considered to affect ice micro- 
physics. Although the parameterization is designed to 
be combined with microphysical schemes that use as- 
sumptions about particle mass, area, and aspect ratio, 
in the case of microphysics parameterizations that do 
not include assumptions about area and/or aspect ratio, 
reasonable assumptions for those quantities can be made 
that are consistent with the assumed mass-dimension 
relations and observations (e.g., Mitchell et al. 1996; 
Baker and Lawson 2006; Auer and Veal 1970; Um and 
McFarquhar 2007, 2009). 

The parameterization described in this paper yields 
the shortwave optical properties of individual ice crys- 
tals that need to be integrated over size distributions. 
Such an integration can be efficiently performed by 
binning the size distribution (Baum et al. 2005b). Since 
the parameterization can yield optical properties for any 
particle size, an unlimited number of bins can be used for 
this numerical integration, although Liu et al. (2012) 
showed that using about 40 bins is generally sufficient. 

To illustrate such an application of the parameteri- 
zation, we calculate fluxes for a uniform cloud layer 
consisting of ice crystals with varying aspect ratios, 
distortion parameters, and projected-area and mass rela- 
tionships. First, we assume mass-dimension and projected 
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area-dimension relationships consistent with aggregates 
of side planes as given by Mitchell et al. (1996); that is, 
a m = 0.0033 gem -6 "’, b m = 2.2, a A = 0.2285 cm 2 '^, and 
b A = 1.88 (such relations are traditionally in cgs units). 
Ice crystal maximum dimension is varied from 1 to 
1000 /um in 1-p.m-wide bins. Note that the power-law 
parameters for aggregates of side planes as given by 
Mitchell et al. (1996) were defined for maximum di- 
mension larger than 600 /im. However, for this illus- 
tration we apply them to all sizes, but limit the mass 
and the ratio of mass and area not to exceed those 
for a solid ice sphere. Furthermore, we consider size- 
independent aspect ratios varying from 0.02 to 50, and 
size-independent distortion parameters varying from 
0 to 0.8. Optical properties are calculated using the 
parameterization for all wavelength bands listed in 
Table 1 and are subsequently integrated over a size 
distribution given by Eq. (20) with /x , p = 1.5 and \ p = 
100 cm 1 (Heymsfield et al. 2002). The resulting bulk 
effective diameter [Eq. (5)] is 84 /im, independent of 
aspect ratio. Note that for some combinations of 
wavelength, maximum dimension, and aspect ratio, the 
particle size parameters may be smaller than the lower 
limit generally accepted for the geometric optics ap- 
proximations applied. However, as discussed above, 
quantification of the resulting errors is currently not 
available. Two-stream calculations of cloud reflectance 
and transmittance are made as described above (see 
also appendix B) for a cloud layer with an optical 
thickness of 4. The cloud layer is embedded between 10 
and 11 km in a tropical atmosphere based on conditions 
observed during the Tropical Warm Pool International 
Cloud Experiment (TWP-ICE) campaign near Darwin, 
Northern Territory, Australia (Fridlind et al. 2012). 
For simplicity, the surface is assumed to be black. The 
results are spectrally integrated to obtain solar up- 
welling fluxes at the top of atmosphere and down- 
welling fluxes at the surface. As shown in Figs. 15a and 15c, 
the variation of radiative fluxes with aspect ratio and 
crystal distortion strongly echo the dependence of the 
asymmetry parameter on these parameters seen in Fig. 9 
(top). Furthermore, the cloud reflectance and trans- 
mittance are almost exactly anticorrelated (not shown) 
leading to a cloud absorptance that is effectively in- 
dependent of aspect ratio and crystal distortion (Fig. 15e). 

Second, to illustrate the additional dependence of 
radiative fluxes and cloud absorptance on projected- 
area and mass relationships, we artificially scale the 
projected area of the ice crystals to obtain bulk effective 
diameters between 25 and 200. [Essentially the same 
results are obtained by scaling masses as apparent from 
Eqs. (5), (21), and (22).] Figures 15b and 15d show the 
dependence of radiative fluxes on effective diameter and 


aspect ratio for smooth crystals (5 = 0). The dependence 
of the fluxes on effective diameter is seen to be weaker 
than the dependence on aspect ratio. Cloud absorp- 
tance, however, is strongly dependent on effective di- 
ameter but only weakly dependent on aspect ratio 
(Fig. 15f), as might be expected. Note that the de- 
pendence of absorptance on effective diameter and as- 
pect ratio resembles the variation of single-scattering 
albedo seen in Fig. 5. Similar dependencies are obtained 
for crystals with greater distortion (S > 0). 

The example seen in Fig. 15 illustrates that ice crystal 
volume, projected area, aspect ratio, and distortion pa- 
rameter that are inconsistent with a model’s micro- 
physics can yield substantial errors in radiative fluxes 
and absorptance (cf. Baran 2012). Our parameterization 
provides for a consistent treatment that avoids such 
errors. Our approach uses the geometric optics ap- 
proximations, which are less applicable for longwave 
radiation calculations. However, longwave radiation is 
mainly determined by absorption, which in turn is de- 
termined by particle volume and projected area, largely 
independent of particle shape (Fu et al. 1999). For long- 
wave flux calculations, optical properties that are largely 
consistent with the shortwave optical properties obtained 
with our parameterization may be calculated by applying 
Lorentz-Mie theory to collections of spheres that con- 
serve ice crystal volume and projected area (Grenfell and 
Warren 1999; Fu et al. 1999; Neshyba et al. 2003), al- 
though such an approach may lead to errors in cloud 
emissivity of order 10% mainly attributable to the over- 
estimation of tunneling effects (Fu et al. 1999; Mitchell 
et al. 2006). The modified anomalous diffraction approx- 
imation (MAD A; Mitchell et al. 2006) partly accounts for 
tunneling effects using shape-dependent tunneling effi- 
ciency factors. 

5. Conclusions 

The parameterization scheme for optical properties of 
individual ice crystals presented in this paper is based on 
the same principles as the well-known scheme devel- 
oped by Fu (1996, 2007), namely, that optical properties 
of complex, aggregated ice crystals can be approximated 
by those of single hexagonal crystals with varying size, 
aspect ratio, and distortion parameter. For example, in 
the case of a bullet rosette, the aspect ratio used in such 
a parameterization is the aspect ratio of the individual 
bullets, while the volume and projected area are the 
quantities determined for the whole rosette. Previous 
schemes by Fu (1996, 2007) and others are restricted to 
preselected wavelength bands, particle geometries, and 
distortion levels. Our scheme further extends these con- 
cepts to provide single-scattering albedos and asymmetry 
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Fig. 15. Illustration of variations of (a),(b) upwelling- and (c),(d) downwelling fluxes and (e),(f) cloud absorptance 
(left) as a function of aspect ratio and distortion and (right) as a function of aspect ratio and effective diameter. The cloud 
layer has an optical thickness of 4 and is embedded in a tropical atmosphere between 10 and 11km. The solar zenith 
angle is 60°. 


parameters for individual ice crystals with any combi- 
nation of volume, projected area, aspect ratio (of par- 
ticle components), and distortion parameter and at any 
wavelength in the shortwave. 


Similar to the scheme developed by Fu (1996, 2007), 
the parameterization of single-scattering albedo and 
asymmetry parameter is based on geometric optics cal- 
culations for single hexagonal crystals with varying size, 
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aspect ratio, and crystal distortion. In the standard geo- 
metric optics implementation used here, the extinction 
cross section is simply twice the projected area of crystals 
in random orientation. The ice crystal single-scattering 
albedo is parameterized in terms of particle aspect ratio 
and an absorption size parameter that is proportional to 
the ratio of the imaginary part of the refractive index to 
the wavelength multiplied by the ratio of particle vol- 
ume to projected area. The parameterization of single- 
scattering albedo uses 26 coefficients. The asymmetry 
parameter of hexagonal ice crystals is parameterized in 
terms of particle aspect ratio, crystal distortion, pro- 
jected area, single-scattering albedo, and wavelength, 
using 62 coefficients in total. 

The scheme is tested for a large variety of hexagonal 
crystals in several wavelength bands from 0.2 to 4 jim, 
revealing absolute errors in both single-scattering 
albedo and asymmetry parameter that are generally 
below 0.015. Furthermore, using a large variety of 
particle aspect ratios, size distributions, and cloud op- 
tical thicknesses, the relative root-mean-square errors 
in cloud reflectance, transmittance, and absorptance 
are shown to be 1.4%, 1.1%, and 3.4%, respectively. 
Note, however, that these errors given do not include 
any possible errors owing to the application of the 
standard geometric optics approximations, as a quan- 
tification of such errors is currently lacking and left for 
future work. 

As we demonstrate, this parameterization scheme is 
flexible enough to obtain ice crystal optical properties 
that are consistent with any assumptions about ice 
crystal mass (equivalent to bulk volume), projected 
area, and aspect ratio, which might be used in modern 
ice microphysics schemes to compute fall speeds and 
depositional growth rates of ice particles. The scheme 
can also be used to estimate optical properties of ice 
crystals of which the mass, projected area, and com- 
ponent aspect ratio are derived from measurements of 
in situ probes or remote sensing (e.g., Mitchell et al. 
1996; Korolev and Isaac 2003; Um and McFarquhar 
2007; Lawson et al. 2010; Mauno et al. 2011; Lindqvist 
et al. 2012; van Diedenhoven et al. 2012a, 2013). Ap- 
propriate distortion parameter values can be derived 
from satellite and aircraft remote sensing measure- 
ments (van Diedenhoven et al. 2012a, b, 2013). 

A Python computer code is provided as supplemen- 
tal material. 


would like to thank Dr. Petri Raisanen for reporting 
a typographical error in Eq. (8) in the early online re- 
lease of this paper. 


APPENDIX A 


Geometrical Relations of Hexagonal Prisms 

Maximum dimension D is equal to 2a hex for plates and 
2«« h ex for columns, where a hex is the side length of the 
hexagon, which is the same as the maximum half-width 
of the hexagonal planes, and a is aspect ratio, determined 
throughout as the ratio of prism height to width. 

Volume V is related to Ohex and a through 

V = 3V3al ex a. (Al) 

Since a hexagonal structure is a convex particle, the 
project area assuming random orientation A p is equal to 
the total surface area divided by 4 (Vouk 1948): 

A p = 3V3al x + l2ai ex a . (A2) 


Effective distance d e , determined through Eq. (4) in 
the main text, is then given by 


, 4 

cl — (It — — . 

e hex v/3 + 4a 


(A3) 


APPENDIX B 


Two-Stream Radiative Transfer 

We apply the two-stream approximation for radiative 
transfer (Coakley and Chylek 1975; Stephens et al. 2001). 
In this method, the reflectance and transmittance of a 
cloud layer with optical thickness t are given by (see 
also Wiscombe and Grams 1976) 


and 


^_ (f/ + l)(L/-l)( e ^o- e -^») 

( u + i)V T A) - ([/ - ^ ’ 


T(H 0 ) = 


4 U 


{U + 1)V t A) - (U - l ) 2 e -rA, 


, (B2) 


Acknowledgments. This material is based upon work 
supported by the NASA ROSES program under Grant 
NNX11AG81G. We are grateful for the contributions 
from three anonymous reviewers. We thank William 
Martin for his help in coding the Python program. We 


where / jl 0 is the cosine of the solar zenith angle and 


U = 


1 — (o + a>( 1 — g) 

1 — a> 


(B3) 
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7 = Vi 1 ~ ")[1 -w + a>{ 1 - g)] . (B4) 

At nonabsorbing wavelengths (&> = 1), this simplifies to 

t(1 - g) 


ROhj) = 


2/x {) + t(1 - g) 


(B5) 


and T = 1 - R. 

The cloud layer can be embedded in a Rayleigh- 
scattering and absorbing atmosphere. Ignoring surface 
reflection and multiple scattering, the monochromatic 
upwelling flux /*' T at the top of atmosphere with a cloud 
layer whose top is located at optical depth t, is given by 


F. = 7? o 

1 t 1 o e 


R T - Itl ndfi, (B6) 


where Fq is the monochromatic solar flux, /jl is the angle 
of outgoing radiation, and R(/x. /jl () ) is a bidirectional 
reflection matrix. To evaluate this integral, we use the 
approximations R(/a, /r 0 ) = R(p- o) and 

f 1 

/xe dfi ~ e ^ T >e t r, (B7) 

. o 

where /z ef( = l/y/3, leading to 

F r = F q R(/x 0 )e _T < (1//i o +1/ ^ff). (B8) 

Similarly, the downwelling diffuse flux F± at the surface 
is given by 


F l = F Q r( / u 0 )e~ T < / ' t o e -( T ,o,-u)A«, (B9) 

where r fc is the optical depth level at the bottom of the 
cloud layer and T tot is the total optical thickness of the 
atmosphere excluding the cloud layer. Note that ex- 
tinction from gaseous absorption and Rayleigh scatter- 
ing within the cloud layer is ignored by Eq. (B9). Cloud 
absorptance a cloud is defined as 

«cioud = [ p o - R (Fo) - n^)]e~ T ^. (BIO) 

For wavelength bands that include high spectrally 
variant absorption lines, the correlated-A: distribution 
method (Lacis and Oinas 1991) is used. 
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